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ABSTRACT 

We present new convergence and stability results for least squares inverse problems in- 
volving systems described by analytic semigroups. The practical importance of these results 
is demonstrated by application to several examples from problems of estimation of material 
parameters in flexible structures using accelerometer data. 
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1. Introduction 

In this presentation we consider abstract inverse problems in a least squares formulation 
for parameter dependent partial differential equations. We are interested in approximation 
ideas which lead to viable computational techniques for such problems. We pursue our 
investigations in the context of the general framework for convergence and stability developed 
by Banks and Ito in [BI] . Motivated by questions related to the use of accelerometer data to 
estimate parameters in flexible structures, we focus on second order (in time) systems with 
sufficient damping so that the system can be modeled by an analytic semigroup. 

We state and prove a new approximation result (a Trotter-Kato type theorem) for ana- 
lytic semigroups. This theorem gives conditions under which a family of approximating semi- 
groups and all its time derivatives converges to a limit semigroup and all its time derivatives, 
respectively. These theoretical results are then stated in terms of simple, readily checked con- 
ditions on the sesquilinear forms defining “stiffness” and “damping” in the abstract second 
order systems. 

We discuss several examples which indicate clearly the practical importance of these new 
convergence results. In two of the examples presented (the damped, cantilevered Euler- 
Bernoulli beam and the “RPL experiment” structure), we can apply the abstract theorem 
to substantially sharpen results already found in the research literature. In a third example 
involving a two dimensional grid structure, we note that the theory can be applied to obtain 
new results for acceleration data convergence in the least squares inverse problems. 

2. Abstract Inverse Problems 

In this section we formulate a class of inverse problems as abstract least squares optimiza- 
tion problems constrained by evolution equations in a Hilbert space and summarize some of 
our previous results for such problems [B], [BI], [BK]. 

We assume we are given a parameter dependent system 

tn i \ *t(0 = i 4(9)' u W + F(t, q ) , 0 < t < T y 

1 > tt( 0) = uo (?) 

for states u{t\ q ) in a Hilbert space H. The parameters q are to be chosen from an admissible 
parameter set Q contained in a metric space (Qi,d). We assume throughout that Q is a 
compact subset of Q\ (a “regularization” assumption). 

We are given observations {z^}, a set of points in the observation space Z, along with an 
observation map C : (7(0, T; i/) — ► Z from the states to the observations. The points z, are 
observations for Cu(U] q), 0 < t; < T, and it is this “data” to which we wish to fit the model 
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by a best choice of the parameter q . Formally, the problem can be stated as 


( v ) 


Find q G Q to minimize over Q the 
functional J(q ) = Y,i \Cu(ti] q) — Z{\z- 


The observation operator C is of fundamental importance to the discussions in this pa- 
per. For parabolic systems (2.1), a typical example for C arises from pointwise evaluation in 
the spatial variables, e.g., Cu{t\q ) = {u(t,Xj\q)} l j=l where Z = R l . For problems involving 
structures such as beams or plates, several examples arise in practice. If one takes measure- 
ments with a laser vibrometer, then one obtains measurements at specific points in space for 
the velocities u t {U t Xj\ q) so that the map C is a composite of time differentiation followed 
by pointwise evaluation. If the measuring devices are accelerometers, one has observations 
for the accelerations u tt (ti, Xj\ q) and thus C is time differentiation (twice) followed by 
pointwise evaluation. 

The inverse problems outlined here are generally infinite dimensional in both the states 
u(t) G H and the parameters q G Q and, moreover, involve unbounded operators related to 
the states ( A(q )) and the observations (C). Thus to develop efficient computational methods, 
one must make finite dimensional approximations for both the state and parameter spaces, 
H and Q, respectively. For the discussions of this paper, we shall restrict our considera- 
tions to approximations of the state space. In some cases (where the parameter sets Q are 
finite dimensional, either naturally or through some a priori parameterization) no approxi- 
mation of the parameter set is necessary. Methods for approximation of Q by sets Q M have 
been fully discussed elsewhere (see [BK]) and these ideas could readily be incorporated in 
our presentation. Since this would add nothing to either the difficulties or their resolution 
that are the focus here, to avoid considerable notational tedium we do not consider such 
approximation ideas in the discussions below. 

Thus, we consider approximate (to be made precise below) state spaces H N C H with 
associated approximate states u N (t ; q) G H N satisfying the approximate systems 

v.»(t) = A N (q)u N (t) + P N F(t, q ) 

( 2 - 2 > u N (0) = P N u 0 {q) 

where A N is an approximation to A and P N is the orthogonal projection of H onto H N . The 
corresponding approximate inverse problems are given by 

N Find q** G Q to minimize over Q the 

' functional J N (q ) = ^^(^i ?) — z %\z- 

Of course, there are a number of reasonable ways in which the approximations may be 
made to arrive at the problems (' V N ). In analyzing different methods and their behavior 
in the context of inverse problems, a number of questions related to parameter convergence 
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arise naturally. For example, given a fixed set of data z = {z<}, do optimal parameters 
of (V N ) converge in some sense to an optimal parameter q for (T 3 )? More generally, one 
might also incorporate the continuous dependence of the estimates on the observations in the 
concept of method stability [B], [BK]. For this, one denotes by q N (z K ) optimal parameters 
obtained from ( V N ) for observations z K and requires that q N {z K ') — » q(z°) in some sense as 
N — » oo and z K — » z° in Z, where q(z°) is an optimal parameter for (V) corresponding to 
observations z°. These issues are carefully discussed in [B], [BK] where it is shown that for 
both parameter convergence and method stability, it suffices to argue that 

For arbitrary sequences {9^} in Q converging 
(2.3) to q € Q, we have Cu N (t\ q N ) — *• Cu(t ; q) as N — » 00 

for each t £ (0, T ). 

Thus, certain fundamental aspects of approximation in inverse problems can be reduced to 
the convergence statement in (2.3) and we shall deal with conditions under which (2.3) can 
be guaranteed in the subsequent discussions of this paper. 

In [BI], Banks and Ito developed a general framework for convergence and stability which 
we summarize here in an extended form to allow treatment of complex valued operators in 
complex state spaces (the sesquilinear forms introduced in [BI] were tacitly assumed to be 
generated by partial differential equations with real coefficients and the theory developed 
there was adequate for such examples). As in [BI], our goal is to isolate properties of the 
system (2.1) and the approximation families to insure that (2.3) holds. We state these 
properties in terms of conditions on a parameter dependent sesquilinear form arising from 
reformulating (2.1) in a weak or variational sense. To this end we consider the system as a 
variational equation holding for all <f> £ V : 

f0 («(<)> <t>) + CT (?)( u (0> 4>) =< ?)> ^ >» 

KlA) t*( 0 ) = 110(g). 

Here V and H are complex Hilbert spaces with V imbedded densely and continuously in 
the pivot space H\ i.e. V *— ► H = H* t — * V *. The sesquilinear form <7(9) : F x F -> C is 
related to the operator A(q) as described below, < . > is the inner product in H, (•, •) is the 
duality product (•, and equation (2.4) is interpreted in the V* sense. The conditions 

on a needed axe: 

(A) Continuity in 9: For 9, 9 £ Q, we have for all <j>, tf> £ V 

- *(§)(& VOI ^ 

(B) V-Coercivity: There exists c\ > 0 and some real Ao such that for 9 £ Q, <f> £ V t we have 

Re p(q)(<f>, <j>) + Ao|^|h > ci|^r. 
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(C) Boundedness in q: There exists c 2 > 0 such that, for q € Q, <f> ) ip € V, we have 


W(q){4>,4 01 < c 2 \4>\ v \Mv- 

Several remarks concerning these conditions are useful. First, we can weaken (A) slightly 
(see [BRR]): From (C) it follows that one can define A(q) in C{V,V*) by cr(q)(<j), ip) = 
(—A(q)^>,xl>)v* > v and then one can replace condition (A) by the requirement that q — ► A{q) 
be continuous from Q to £(V, V*). 

Moreover, for a : V x V —* C defined using real coefficients, i.e. real scalars, it suffices 
to replace (B) by the condition: cr(q)((f>,<f>) + A 0 |^|jj > cx|^|fr for all real valued functions ^ 
in V (this is condition (B) as stated in [B I] ) . To see this, we note that cr(q)(<f) + iip, $ + iip) 
= <Kg)(& <t>) + VO + » W9)(^> 4>) - VO] and thus Re a{q){<t> + V + iVO = 

cr(q)(<f>,<f>) + cr^XVs^). Therefore, a (q)(<j>,<f>) + ^o\4>^h — c il^lv f° r a ^ rea l valued <f> mV 
implies immediately that condition (B) above holds. 

Under conditions (B) and (C), a defines operators A(q) such that cr(q)(<f> } i{)) = 
< —A(q)<}>,Tp > for <f> (E dom (A(q)),‘ip (E V with dom(A.(g)) dense in V (e.g., see [S]). More- 
over, A(q) is a sectorial operator with (A I — A(g))dom(A(g)) = H for all A with ReA > A 0 ; 
indeed Rx(A(q)) = (A I — A(g)) -1 exists as a bounded operator on H for all A in the com- 
plement of a sector with A 0 as a vertex. The operator A(g) generates an analytic semigroup 
T(t\q ) which can be used in defining mild solutions of the system (2.4) or (2.1); i.e., 


(2.5) 


u(t\q) = T(t\ q)u 0 (q) + f T(t - s; q)F(s, q)da. 

Jo 


Galerkin approximations (e.g., systems (2.2)) for (2.5) can be developed in the context of 
sesquilinear forms satisfying conditions (B) and (C). Let H N be a family of finite dimensional 
subspaces of H satisfying H N C V and the condition 


(Cl) 


For each <f> € V, there exists (j> N € H N such that 
\<f> — F\ v — > 0 as N — ► co. 


(We note that many popular approximation schemes - e.g., linear splines and more generally, 
many finite element schemes - satisfy condition (Cl).) To define A N (q) : H N — > H N , we 
restrict a(q) to H N x H N and denote this restriction by cr N (q)- Then a N satisfies conditions 
(A), (B), (C) on H N x H N which implies existence of a bounded linear operator A N (q) 
on H n satisfying <J N (q)(<j>, VO =< ~-A N (q)<f> ) ip > for all <f>,ip £ H N . From (B) and (C) it 
follows that the A^(g) are uniformly dissipative ^sectorial) and generate analytic semigroups 
T"((; q) on H N . Solutions of the systems (2.2) are then given by 

(2.6) u K (t 1? ) = T N ((; ,)?"%(,)+ f T !, (t-s;q)P' , F(s,q)da. 

JO 
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It can then be established that under conditions (A), (B), (C) and (Cl), we have 
u N (t\q N ) — > u(t\q ) in V for arbitrary sequences { 9 ^} with q N — * q in Q. Observe that 
this immediately yields (2.3) if the operator C possesses certain boundedness properties. 
The focus of this note involves cases where the operator C does not have such boundedness. 

The arguments for the convergence u N (t; q N ) — > u(t] 9) for the situation here (complex 
Hilbert spaces V y H and complex valued sesquilinear form a) are essentially the same as 
those in [BI], Slight changes in the arguments for Theorem 2.2 of [BI] are necessary to treat 
the case of complex valued inner products and sesquilinear forms (essentially one need only 
use Re cr(-, •) and Re < •, • > in some of the inequalities). In these arguments one does not 
use directly the analyticity properties of the semigroups T(t]q),T N (t]q). Rather one relies 
heavily on resolvent estimates of Tanabe [T] depending on the V-coercivity of a along with a 
resolvent convergence form of the Trotter-Kato approximation theorem of linear semigroup 
theory (see [BI] for details). 

Our interest here is mainly in second order systems of the form 

u(t) + B(q)ii(t) + A(q)u(t) = f(t,q ) 

(2.7) u(0) = u 0 

u(0) = uo 

where A(q) is a generalized stiffness operator and B{q) is a generalized damping operator. 
Again we consider this equation in a weak or variational sense defined via parameter depen- 
dent sesquilinear forms in a complex Hilbert space V t — » H = H* > V*. We are given a 
stiffness sesquilinear form cri(q) : V xV — > C that is symmetric and satisfies the boundedness 
condition (C). Then there exists A(q) G C{V y V*) such that o\ (g)(^,V0 = (A(q)<f>ti>)v,v- 
We also have a damping sesquilinear form £ 7 - 2 ( 9 ) : V X V — > C which satisfies (C) so that 
there exists B(q) G C(V,V*} with 0- 2 (g)(^,VO = {B(q)<f>,'tp)v*,v- We then reformulate the 
system (2.7) and seek solutions u(t) G V satisfying for all $ G V 

, 0 (*(<)> 4>) + ° 2 (q)(u(t), <f>) + o-i (g)M*), <t>) =< f(t, q), <t> > 

' ' ' u( 0 ) = u 0 , ii( 0 ) = v 0 . 

As is standard practice, we rewrite this in first order vector form on = V x H and 
V = V x V in the coordinates (it,u). To this end, define a(q) : V x V — * C by 

(2.9) <Kg)(K i>), (<f>, VO) = - < v, <f> >v +< 7 i( 9 )(u, VO + a 2 (q)(v, VO 


so that (2.8) may be rewritten as 


W).X) + ^(?)(4).X) =< B(t,q),x >H 
u;(0) = (u 0 ,v 0 ) 
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for w(t ) = (u(t),u(t)) and x = (^,V0 i n V with F(t,q ) = (0 ,f(t,q)). Or, if as in the usual 
practice, we abuse notation and do not distinguish between row and column vectors, we may 
write this in equivalent operator form 

w(t) = A(q)w(t) + F(t,q) 
u;(0) = (u 0 ,v 0 ) 

where cr{q){x, 0 = with 

,, v F o i 

A{,) - [ -4( ? ) -B(q) 

For our treatment here (as for that in [BI]) we assume that a* (q) satisfies conditions 
(A), (B), and (C). We further assume that a 2 (q) satisfies conditions (A) and (C). Then the 
strength of the coercivity assumption on < 72 ( 9 ) determines the properties of the semigroup 
generated by A(q). For example, if <7 2 (<?) satisfies (B), then *4(g) is strongly V-coercive and 
generates an analytic semigroup on 7i. (Actually, A(q) is V,-coercive with V, = V q x V q 
where is V taken with the equivalent inner product <7i(g)(v) — see [BI].) In [BI], it is 
assumed only that is W-semi coercive: There exists b > 0 such that for all <f> € V we have Re 

t7 (9)(^> 4>) b\4>\)i- I 11 this case one can argue only that A(q) generates a strongly continuous 

semigroup on H. The theory for second order systems with this weak damping is developed 
in [BI] in order to treat several forms of damping (spatial hysteresis, time hysteresis, bending 
rate damping) which are of physical interest and yet do not satisfy the strong V — coercivity 
assumption. The convergence theory obtained yields that u N (t\ q N ) — > u{i\q) in V norm, 
u N {t] q N ) — ► ii(f ; q) in H norm whenever q N — » q. 

As opposed to [BI], we wish to consider in this paper the case where a 2 (q) does satisfy 
the strong V — coercivity condition (B) and hence A(q) is the infinitesimal generator of an 
analytic semigroup. We obtain immediately that cr(q) of (2.9) is V— coercive and the first 
order theory outlined above can be applied directly to the system (2.10). This yields the 
convergence statement in V - i.e., u N {t\q N ) — ► u(t]q) in V , u N (t]q N ) — » u(t',q) in V. A 
natural question arises as to whether we can use the analyticity of the semigroups in this case 
to obtain stronger results. The next section is devoted to results that yield an affirmative 
answer to this question. 

3. Analytic Semigroups and Approximation 

We first give a general approximation theorem for analytic semigroups that is a general- 
ization of the well-known Trotter-Kato theorem [BK]. 

Theorem 3.1. Suppose we have complex Hilbert spaces X and X N , N = 1,2..., with 
X N C X. Let P N : X — > X N denote the orthogonal projection of X onto X N satisfying 



6 


p N J strongly. Suppose that A N and A are the infinitesimal generators of analytic 
semigroups S N (t ) and S(t) on X N and X respectively that satisfy the following: 

There exists a region £ = £$ = {A £ C : |ar</(A — Ao)| < f + $}, where 6 > 0, 

such that £ U{A 0 } C p(A) D^ =1 p(A N ) and 

(i) there exists a constant M independent of N such that 

for all A £ £ and N = 1,2,...; 

(ii) for some A £ £ and each x £ X we have R\(A N )P N x — > Rx{A)x. 

Then we have 

(Hi) For each x € X, S N (t)P N x — * 5(t)x uniformly in t on compact 
subintervals of [0, oo); 


(iv) For each x £ X and integer k > 1, (A N ) k S N (t)P N x — > A k S(t)x 
uniformly in t on compact subintervals of (0,oo). 

Proof. The statement of (iii) under the given conditions is just a variant of the well-known 
Trotter-Kato theorem and follows immediately from Theorem II.1.14 of [BK]. To argue result 
(iv), we first note the convergence in (ii) for some A £ £ implies that the convergence holds 
for all A £ £ . In light of the resolvent bounds of (i), this follows from the identity for 

A £ £ 

r^A n )P n - RM) = [/ + (A - p)R ti (A N )P N )[Rx(A N )P N - Rx(A)][I + (A - fi)R^A)) 

which is readily established using the standard resolvent identity R\(A) — R M (A) = (p — 
A )R x (A)R m (A). 

The analyticity of the semigroups S N (t ) and S(t) allow us to write for t > 0 (e.g., see 

[P]) 

A t S(t) = -.J r X^A)dX 

and 

(A N ) k S N (t)P N = -L f \ k e Xt Rx(A N )P N d\. 

L'KX J P 

Here T is a positively oriented contour through Ao lying in £ with arg(A — A 0 ) = iv for 
A ^ A 0 where v is fixed in (^-f 8— e, ^+£)- The desired convergence results follow immediately 
from the inequality 

| (A N ) k S N (t)P N x - A*S(f)x| 

<f [ |A|*|e i *| | R x (A n )P h z - JJ;>(A)z|<iA 
27 t Jr 
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using the resolvent convergence of (ii) by noting that the integrand in this integral is domi- 
nated by an integrable function. 

We now return to the second order systems of Section 2 - see (2.8), (2.9), (2.10) - and 
use Theorem 3.1 to obtain our main convergence results. As before we take Ti —V x H and 

V = V x V. Let Ti N = H N x H N and P N be the orthogonal projection of Ti onto Ti N . 

Theorem 3.2 Let &i(q) and cr 2 (q) in (2.8) satisfy conditions (A), (B) and (C) and let 
H N C V satisfy condition (Cl). Let {g^} be arbitrary in Q with q N — ► q. Then we have 

(i) The sesquilinear form cr{q) given by (2.9) satisfies conditions (A), (B), (C) in the norms of 

V andTt and the operator A(q) defined via a (g)(x,£) =< —A(q)x,t >n for x £ dom{A(q)) 
is the infinitesimal generator of an analytic semigroup T(f;g) on Tt. 

(ii) Let A N (q) denote the operator obtained by restricting a(g) to TL N x Ti N and let T N (t\ q) 
denote the corresponding analytic semigroups on Tt N . Then we have 

(a) For each x € Tt, T N {t\ q N )P N X q)x in Tt uniformly in t 

on compact subintervals of [0, oo); 

(b) For each x € Ti and positive integer k, A N (q N ) k T N {t\q N )P N x —> A{q) k T{t\q)x 
in Ti uniformly in t on compact subintervals o/(0, oo). 


The proof of this theorem is now rather straightforward. We first note that it is readily 
shown with routine calculations that oi,< 7 2 satisfying (A), (B), (C) in the V and H norms 
imply that cr satisfies (A), (B), (C) with the norms of V = V x V and Ti = V x H. Result 
(i) of the theorem thus is established. 

For result (ii), we apply Theorem 3.1 with X = Tt = V x H and X N = Ti N = H N x H N 
along with arguments of [BI]. Restricting cr{q N ) to Ti N x Tt N to obtain the operators A N = 
A N (q N ), we have that condition (B) is satisfied in Tt N with the constants uniform in N - 
i.e., the uniform sector condition and the uniform resolvent bounds in (i) of Theorem 3.1 are 
readily seen to hold. For the resolvent convergence of (ii) we refer to Theorem 2.2 of [BI], 
noting that we have all the hypotheses of that theorem holding here in the sense of the V 
and Tt norms. Hence the same arguments (modified slightly as mentioned in Section 2 above 
to account for complex Hilbert spaces V and Tt) given for Theorem 2.2 of [BI] can be used 
here to establish (ii) of Theorem 3.1. The convergence statements (a) and (b) of Theorem 
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3.2 then follow directly from (iii) and (iv) of Theorem 3.1. 

We are now in a position to use the results of Theorem 3.2 to address the question of 
stronger convergence results for u N (t\ q N ) to u(t\ q) raised in connection with the observation 
operator C in Section 2. We first consider the case of homogeneous systems (i.e., / = 0 
in (2.8)) in which case w N {t-,q N ) = {u N {t]q N ) } u N {t]q N )) = T N (t-,q N )P N X and w{t-,q) = 
{u(t-,q),u{t\q)) = T(t;g)x • For x <*om (.4(g)) we have u;^(t) = A N {q N )T N {t\q N )P N x, 
w"(t) = A N {q N ) 2 r N {t-,q N )P N X , etc., while w(t) = A{q)T(t] q) X , w{t) = A(q) 2 T (t) q) X , 
etc. Thus, from statement (b) of Theorem 3.2 we find w N (t]q N ) -* w(t\q) in H, which 
yields u£(i; q N ) -»■ u H (t] q) in the V norm, uniformly in t on compact subintervals of (0,oo). 

Indeed, for all k = 1, 2 . . . , we obtain £u N (t; q N ) -» j£u(i; g) in V. 

Since, as we shall see in the next section, many interesting examples involve V C H (fl) 
with Cl C R l or fi C R 2 , and since in this case V imbeds continuously in C{Cl) t the above 
results guarantee pointwise (in t and x) convergence of u N (q N ) and all its time deviatives 
to tt(g) and its time deviatives, respectively. Thus, observation operators related to laser vi- 
brometers (ut) and accelerometers (u tt ) are included in the convergence and method stability 
theory for output least squares inverse problems. 

Similar results are available for the nonhomogeneous equation (2.8) (or (2.10)) if one 
places appropriate regularity conditions on /. To obtain these, one considers the representa- 
tion theorems (mild solutions) 

w N {t) q N ) = q N )P N w 0 + Jo* T N {t - s; q N )P N F{a, q)ds 

w[t ; g) = T (t; g)u;o + So T{t - s] q)F(s , q)ds 

and uses regularity conditions on mild solutions in 7i = V x H. For example, if f(‘,q) is 
in C7 1 ([0,T],if), one can differentiate once, if /(-,g) is in C 2 ([0,T],H), one can differentiate 
twice, etc. (see [P]). Indeed, since we have analytic semigroups, we can weaken the conditions 
on / needed for this procedure by considering special regularity theorems for mild solutions 
of Cauchy initial value problems in the case of analyticity (e.g., see Chapter 4.3 of [P]). For 
example, it suffices to have /(-, g) G L l (Q,T]H) and locally Holder continuous in order to 
differentiate w ^ and w once in the above representations. 

4. Examples 

In this section we present briefly several examples to which the above theory can be 
readily applied. In some cases this sharpens the convergence and stability results currently 
available in the literature; in other cases it provides new results for the associated inverse 
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problems. 


Example 4.1. We consider a cantilevered Euler-Bernoulli beam with Kelvin- Voigt damping 
(i.e., stress proportional to a linear combination of the strain and the strain rate). The beam 
is assumed fixed at x = 0 and free at x = l, with the transverse displacement at time t and 
position x given by u(t, x). Typical observations consist of acceleration u u (t, Xj) or velocity 
u t(^> several locations Xj. Balance of forces and moments yield the following system 

(we assume the linear mass density p is normalized to unity): 

d 2 u d 2 , d 2 u d^u 

W + + CDl axm^ = 0<x<t > 


(4.1) 


Oil 

u(t,o)= ^(f,o) = o, 


r d 2 u d z u\ 

^ EI dx » + CDl dx^ x=l = °’ 





(Pu 

dx 2 dt 


}).=< = 0. 


The parameters to be estimated in typical examples (see [BWIC], [BFWIC], [Bln]) are the 


stiffness and damping coefficients, El and cq/, respectively. For a parameter set we choose 
Q compact in 


Q v = {q = ( EI,c d I ) : q 6 L°°(0,£) x EI(x) >v> 0 ,c D I(x) > v}. 


For the state spaces H and V we take H = H°(0,£), V = Hl(0,£) = {<f> E H 2 (0 ,£ ) : ^(0) = 
^(0) = 0}. Then in writing the system (4.1) in the form (2.8) we take 


<r 1 (q)(<f>,iP)=<EID 2 <t>,D 2 i> >o 
*2{qMTp)=<c D ID 2 <l>,D 2 iP > 0 


where D 2 = and < •,• >o is the inner product in H°[0, £) = L 2 (0,£). For real valued <f> 
in V we have at once 


4>) > v\D 2 <j ) | 0 > Cx\<j>\\ 


which, as we have noted, implies 


Re ai(q)(<f>, <f>) > c x \<j>\%- 


for complex valued <f> in V. It is equally trivial to argue that conditions (A) and (C) of 
Section 2 hold for this c^q), o^q). For Galerkin schemes satisfying condition (Cl) we thus 
have Theorem 3.2 applicable, and hence parameter convergence and method stability hold 
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iHM III 



for the least squares problems (T 3 ) and (^P N ) whenever one formulates these with pointwise 
observations of either velocity or acceleration. 


Example 4.2 For this example we return to the so-called “RPL experiment” discussed in 
some detail in [BGRW]. The focus of our attention is a cantilevered Euler-Bernoulli beam 
with a flexible gas hose and thruster nozzle attached to the free end as depicted in Figure 
2.2 of [BGRW]. The structure is modeled as a uniform cantilevered beam with Kelvin- Voigt 
internal damping and tip mass with a mass-spring-dashpot assembly attached at the tip. 
Along with the usual damped Euler-Bernoulli beam equation 


(4.2) 


d 2 U , r.id'u , T d‘u 

+ a? + CDl a7ai 


= 0, 0 < X < l, 


we have the force balance equation at the tip 
(4.3) 


d 2 u d*u 

~ CoI d^di 


ti 

" EI W ] ~‘ 


= <*( y(t) - + W»(0 - Mt, 0) + /( 0 


and the hose assembly state equation 


(4.4) 


3il 

m H y(t) + cn(y(0 - “ «(*.*)) = °* 


Here p,EI } and cpl are the usual beam parameters whereas m? represents the tip mass, 
m H is the hose mass, k H is the hose stiffness (the “spring” or restoring force constant), c H 
is the hose damping coefficient and /(<) represents an externally applied force at the tip 
(firing of the tip mounted thrusters). Boundary conditions for the coupled state equations 
(4.2) and (4.4) include the tip force balance equation (4.3), the tip moment balance equation 
(assuming that the hose assembly has negligible rotational inertia) 


(4.5) 


d 2 u (Pu . 

I BI a? + CDl B^m ] “‘ = °> 


and the zero displacement, zero slope conditions at the fixed end x — 0 


(4.6) 


Q 

u(f,0) = ^(t,o) = °* 


The structure is assumed initially at rest so that initial conditions are given by 

8il 

u(0,x) = -^-(0, as) = 0 

y(o ) = y(0) = 0. 
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The parameters to be estimated using accelerometer observations (see [BGRW]) include 
q = (ray, El, cpI, mu, ch, kn), which is to be chosen from a compact subset Q C R+ ■ 

To write the system (4.2)- (4.6) in weak or variational form, we use the state variable 
u(t) = (y(t),u(t,£),u(t, •)) in the state space H = R 2 x H°( 0,1) with inner product 

< (C> v, 4> ), (A, n, i>) >h=~& + rni+ < <f>, i; > 0 . 

For the space V we choose V = G H : <f> G H 2 (0,£), ^(0) = D<f>(0 ) = 0, 77 = <}>(£)} 

with inner product 

< (c. m *), (a, m,*) >v= (c- m )( a - m)+ < >» . 

The stiffness and damping sesquilinear forms are given for $> — ((, <j>{£), <j>), = (X,‘ip(£),i>) 

in V by 

= k H (C - m ){ A - VW) + El < D 2 4,D 2 i > 0 
0 - 2 (?)(<£, = c/r(C - ^(i))(A - V»(£)) + c D I < D 2 <f>, D 2 ip > 0 . 

We also need an operator M(q) G £(77) given by 

•M(g)(C,7?, = (,'m H C,m T rj,p<f>). 


This operator can be extended to C{V*) in an obvious manner. Then the system (4.2)-(4.6) 
can be written in variational form for the state u(t) G V to satisfy for all ^ G V 


(M(q)utt(t)>4>) + ffa(9)(«t(0> i) + 01 {q)(u(t),j>) =< F(t),J>> 

(4.7) 

*(0) = u*(0) = 0, 

where F(t ) = (0,/(t),0). Since, for mH,rriT,p positive, the operator M(q) is invertible, 
equation (4.7) is obviously equivalent to an equation of the form (2.8). Thus the theory of 
Sections 2 and 3 is applicable if <j\ and Vi satisfy the requisite hypotheses. If Q is bounded 
below in 72®, it is readily seen that both U\ and 03 satisfy conditions (A), (B), and (C). For 
example, we see immediately for <f> = ((, <f>(£), <f>) G V 


*>(»)(#, i) = chU - M) 2 + od1\d^\1 

> d{(( - m f + \D 2 t\l} = C t \i\ 2 y. 

Similar arguments hold for ai(q ). 

Application of the theory in Sections 2 and 3 substantially sharpens the results given in 
[BGRW]. In that paper the main results (see Lemma 3.1) yield 

(4.8) / |ujf(t; q N ) - Utt{t] q)\ H dt -* 0 

Jo 
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a8 jy — » oo. The arguments are rather tedious and require the assumption ii(q) € H 2 ( 0 , T\ V) 
on the limit function. Note also that (4.8) would require continuous time acceleration ob- 
servations be used in the least squares criterion. 

In contrast, the theory of this paper yields (since H 2 embeds compactly in C ) u N (t, x; q N ) -> 
u(t,x-,q), u N (t,x\q N ) — > u (t,x;q) uniformly in x G [0,£] for each t € [0,T] as well as 
u%(t\q N ) -> Uu{t]q N ) in V for each t, which permits sampled time acceleration observa- 
tions. The arguments are simple (given the theory developed above) and do not require the 
a priori smoothness assumption on u(q). 


Example 4.3 As a final example, we briefly describe the models for two dimensional grid 
structures (“plates with holes”) developed and investigated in [R] and [BR]. We use Love- 
Kirchoff plate theory with Kelvin- Voigt damping. The “plate” is rectangular, e.g., on (x,y) 6 
[0,€ a ] x [0,4], with rectangular holes periodically placed to produce a thin planar grid. It is 
assumed to be hanging vertically, clamped at the top, with the other three edges free. The 
basic equation for transverse displacements w(t,x,y ) is given by 


d 2 w d 2 M x n d 2 M xy d 2 M v _ 
^H 2 + ~d^~ + 2 dxdy + dy 2 


where h is the thickness of the plate, the bending moments Af®, M v about the x and y axes, 
respectively, are given by 


M x - 


El r d 2 w d 2 w , c D I . d 3 w 


1 — v 2 1 dx 2 


{ a? + " + 


dx 2 dt 


+ v 


El r d 2 w d 2 w . c D I { &w 
MV ~ 1 - i/ 2 ^ ~dy 2 +V dx 2 ^ + 1 -v 2 ^dy 2 dt + U 


d 3 w 

dy 2 dt 

d 3 w 

dx 2 dt 


} 

}. 


and the twisting moment is given by 


EI f d 2 w d 2 w c D I . d 3 w _ d 3 w 

^ 1 — v 2 ^ dxdy dxdy 1 — v 2 dxdydt dxdydt 


Here v is Poisson’s ratio. For the boundary conditions at the top (clamped) we have the 
essential boundary conditions w = ^ = 0 (the top corresponds to the x— axis). The plate 
is free on the other outer edges and on the edges of the holes, where the natural boundary 
conditions of zero moment and zero shear are required. For example, on a free edge parallel 
to the y — axis this results in the conditions 


M* = 0, 


dM x 2 dM xv 
dx dy 


As shown in [R], [BR], one can readily define the corresponding stiffness and damping 
sesquilinear forms Cx^Oi on V x V , where V — G H 2 (il) : <f> = ^ = 0 along y = 0}, fl 
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is [0,4] x [0,4] l e8S the holes, and argue the needed V —coerciveness along with conditions 
(A) and (C). In this case, the state space is H = The theory of Sections 3 and 4 can 

thus be shown to hold for this example. 
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